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Abstract -We report on the effect of substrate temperature (T) on both local structure and 
long-wavelengtli fluctuations of polycrystalline CdTe thin films deposited on Si(001). A strong 
T-dependent mound evolution is observed and explained in terms of the energy barrier to inter¬ 
grain diffusion at grain boundaries, as corroborated by Monte Carlo simulations. This leads 
to transitions from uncorrelated growth to a crossover from random-to-correlated growth and 
transient anomalous scaling as T increases. Due to these finite-time effects, we were not able to 
determine the universality class of the system through the critical exponents. Nevertheless, we 
demonstrate that this can be circumvented by analyzing height, roughness and maximal height 
distributions, which allow us to prove that CdTe grows asymptotically according to the Kardar- 
Parisi-Zhang (KPZ) equation in a broad range of T. More important, one finds positive (negative) 
velocity excess in the growth at low (high) T, indicating that it is possible to control the KPZ 
non-linearity by adjusting the temperature. 


Thin films are the basis of the optoelectronic indus¬ 
try. Commonly, patterned/mounded interfaces are ob¬ 
served due to growth instabilities [T] or polycrystallinity, 
where a complex growth involving intra- and inter-grain 
dynamics arise. It is well-known that size, texture and 
spatial distribution of these structures affect several thin- 
film properties that are crucial for applications in solar 
cells [2], spintronic devices [3], contact technology [3] and 
many others. 

At a coarse-grained level, the evolution of thin films 
and other growing interfaces is also a subject of broad 
interest, since they exhibit scaling invariance and univer¬ 
sality [SM]- For instance, the kinetic roughening of flame 
fronts [7], turbulent phases in liquid crystals [ 8 ], colloidal 
particles deposited at the edges of evaporating drops [9], 
silica ina, CdTe m and oligomer films m have been 
shown to belong to the celebrated Kardar-Parisi-Zhang 
(KPZ) pj] universality class (UC). Thus, despite their dis¬ 
tinct microscopic nature, the interface of all these systems 
evolve asymptotically according to the KPZ equation m- 

9h d*’^ = iA7 2 /i+ ^(Vfr) 2 +r)(x.,t), (1) 


where h(x, t ) is the height at substrate position x and time 
t, and v, A and rj account, respectively, for surface tension, 
interfacial velocity excess and white noise. 

Since the design of thin films with specific properties 
requires control of their growth, a natural question raises 
up: How do the growth parameters affect the local and 
long-wavelength dynamics of the system? In fact, the ef¬ 
fects of parameters such as substrate temperature El, 
molecular flux [13] and electric potential m have already 
been studied. However, the (roughness) dynamic scaling 
analysis (DSA) performed there did not lead to any con¬ 
clusion about the UC of those systems, possibly due to 
corrections to scaling/transient effects. 

In this Letter, we answer that question regarding the 
effect of deposition temperature (T) on the growth of 
CdTe thin films, a very important material for the fabri¬ 
cation of detectors (of 7 - and X-rays) [14] , solar cells [17] , 
ultra-fast optical sensors [IS] and others [T9] . The struc¬ 
ture and morphology of films, grown by Hot Wall Epitaxy 
(HWE) at different T, have been characterized by X-ray 
diffraction (XRD) and atomic force microscopy (AFM). 
A complex mound evolution is observed and has been ex¬ 
plained in terms of diffusion barriers at grain boundaries, 
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Fig. 1: 10/im x 10/rm AFM images for films grown at T = 200°C [(a) and (b)] and T = 300°C [(c) and (d)] for 120 (left) 
and 240 min (right). Cross sections of typical mounds for T = 200 and 300 °C are shown in (e) and (f), respectively, (g) XRD 
results for T = 250 °C and deposition time t = 240min. The inset shows the coefficient of preferential growth in [111] direction 
as a function of growth time. 


as confirmed by Monte Carlo simulations. Although tran¬ 
sient/crossover effects in the DSA have also been found 
here, through the study of several distributions, we prove 
that long-wavelength dynamics of the CdTe surfaces be¬ 
long to the KPZ class in a broad range of T. The In¬ 
dependence of the parameters in Eq. [T] is also unveiled. 

Experimental methods. — CdTe (5N) was evapo¬ 
rated on Si(001) substrates by HWE, a well-controlled 
and highly reproducible growth technique [2DJ1H]- The 
experimental HWE setup is described in detail in Ref. 
m • Substrate cleaning and growth conditions are the 
same reported in m, however, temperatures T = 150, 
200 and 300 °C are also considered here. The HWE tech¬ 
nique has been chosen because the growth apparatus is 
much simpler and the running costs are lower than those 
for Molecular Beam Epitaxy (MBE). The growth rate was 
determined ex-situ measuring the sample thickness with a 
ContourGT-K (Bruker) optical profiler. Surface topogra¬ 
phies were measured in air using a Ntegra Prima (NT- 
MDT) SPM in contact mode with Si tips. Images of 
10/im x lOfxm (1024 x 1024 pixels) were carried out for 
3-10 distinct regions near the film center. Crystallinity 
and texture features were investigated by X-ray diffrac¬ 
tion (XRD) using a Bruker D8-Discover diffractometer. 

Local dynamic. — Figures [lji-d show typical AFM 
images for CdTe films grown at T = 200 and 300 °C 
for 120 and 240 min. The grained/mounded morphology 
expected for polycrystalline surfaces is observed, where 
for T = 200 °C conical grains with well-defined bound¬ 
aries dominate the surface at short times (Fig. [ft). For 
the largest time available (Fig. [ft) one sees the pres¬ 
ence of some coalesced/packed graintQ carrying a multi- 

1 The large structures (mounds) can be formed by both coales¬ 
cence of grains with the same crystallographic orientation and pack¬ 
ing of grains with different orientations. 


peaked form. There is a plenty of these structures at sur¬ 
faces grown at T = 300 °C since short growth times (Fig 
m , and they keep growing to give place to large mounds 
separated by deep valleys (Fig ftl). These features are 
highlighted in Figs, [ft and [ft, where cross sections of 
characteristic superficial structures are shown. The poly¬ 
crystalline environment is confirmed by the appearance 
of several peaks in the 9 — 29 XRD spectra. It is pre¬ 
sented in Fig. [ft for T = 250 °C, and similar spectra are 
found for all investigated T. Additionally, a strong (111) 
T-independent texture is revealed (see inset of Fig. [ft), 
pointing out that (111) grains grow faster than the oth¬ 
ers and, upon coalescence/packing, cover the neighboring 
non-(lll) ones. Previous studies suggest that this CdTe 
texture is also independent of the substrate [22] . 

Figures [ft and ft show the local roughness [wi oc (l,t) - 
defined as the rms height fluctuation inside a box of lat¬ 
eral size l] versus l for T = 200 and 300 °C, respectively. 
Solid lines indicate linear fits used to extract the local ex¬ 
ponent ot \, defined by wi oc ~ l ai [23]. One notices that 
this exponent is measured for l < 0.1/rm, so that it charac¬ 
terizes the intra-mound morphology. As demonstrated in 
Ref. [23], 04 decreases, as sharper are the mound shapes, 
from ou ss 0.90 (for rounded mounds) down to a\ « 0.50. 
Therefore, the exponents depicted in the inset of Fig. [ft 
indicate that, for a given T, the top of mounds becomes 
more rounded as time evolves. The same is seen when 
the time is fixed and T increases, as corroborated in Figs. 
[TJi-f. Unfortunately, a crossover to the truly roughness ex¬ 
ponent [23] is not observed in the local roughness (Figs, 
ft-b) . 

One may notice in Fig. [ft, at short-length scales (/ < 
0.1/irn), wio C increasing in time for t ^ 120 min. This is 
the hallmark of anomalous scaling [ft] , but this “anomaly” 
is transient, since at large t one sees Wi oc decreasing in 
time, leading to the standard Family-Vicsek scaling [25] . 
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Fig. 2: Local roughness wi oc versus box size l for (a) T = 
200 and (b) 300 °C. Solid lines indicate linear fits used to 
extract the local exponent ot\. In (a) and (b) insets show qi 
and {(VI) 2 ) against time, respectively, for T = 150 °C (blue 
triangles), 200 °C (green diamonds) and 300 °C (red circles). 

A similar behavior was found in CdTe films grown at T = 
250 °C m- The time evolution of Wi oc at short-length 
scales is directly related to the spatially averaged squared 
local slopes ((V/i) 2 ) at interface (see the inset of Fig. [5b). 

The origin of these results can be understood as fol¬ 
lows. Initially, CdTe grains evolve in the Volmer-Weber 
growth mode [5B], and as higher is T larger are their 
widths [D (Fig. EJ. Moreover, one found that the grain 
aspect ratio (r = height/width) also increases with T at 
short times, possibly due to an unbalance between up- and 
downward diffusion at grain edges. As the initially isolated 
grains enlarge laterally, they collide forming grain bound¬ 
aries (GBs), where defects are formed. These defects give 
rise to an additional energy barrier Egb to diffuse toward 
these sites m, as also suggested recently in the growth 
of CdTe/CdS films [55]. At low T, a small number of 
molecules overcomes this barrier and most of them ag¬ 
gregates inside the grain where they have arrived. This 
compels the grain height to increase faster than its width, 
leading r and ((V/i) 2 ) to increase (Figs. and [ 2 ) 3 ) . As 
time evolves, aggregations at GBs induces a relaxation, 
which diminishes the number of superficial defects [29] 
and, consequently, the inter-grain diffusion becomes more 
active, since the Egb barrier disappear in those relaxed 
places. Thence, the coalescence/packing of grains becomes 
more operative and small grains give place to larger width 
structures, so that r and ((V/i) 2 ) start to decrease, ex¬ 
actly as observed for T = 200 °C (Fig. [5b). For higher T, 
where the surface diffusion is more active, the relaxation 
process happens earlier as well as the decreasing in r and 
((V/i) 2 ) (Figs. QJ, [2jo and Ref. ED- 

The reliability of the above reasoning is illustrated in a 
very simplified one-dimensional atomistic growth model. 
Since our interest is the coalescence process, the growth 
starts on a periodic array of pyramidal grains with the 
same width £ and height H , for simplicity. A randomly 
deposited particle diffuses at surface until reaches a site 
i satisfying the constraint | hi — hi± i| < 1, where it per- 



Fig. 3: (a) Height profiles for T = 150 (top) and 200 °C (bot¬ 
tom) for t = 10, 100 and 1000, and shifted by k = 10, 80 
and 960, respectively. The vertical lines indicate the location 
of the initial GBs. (b) Average squared local slope {(Vfo) 2 ) 
versus time. 

manently aggregates. Therefore, inside the grain, aggre¬ 
gation follows the conservative RSOS (restricted solid-on- 
solid) rule [3D]. However, at the GBs there is an energy 
barrier Egb-, so that a particle diffuses toward them with 
probability Pjj = e~ EaB / kBT . Once a particle aggregates 
at a given GB i, the barrier Egb at i becomes null with 
probability P R = e~ ER / kBT , in order to mimic the relax¬ 
ation process. Figure [3b shows typical surface evolutions 
for T = 150 and 200 °C, with E G b = 0.10 eV, E R = 0.30 
eV, C = 64 and H = 8, 16 and 24 for T = 150, 200 
and 300 °C, respectively. For T = 150 °C one observes 
grains with almost fixed widths and increasing heights. A 
similar behavior is found at short times for T = 200 °C 
but, for long t, large mounds (formed by coalesced grains) 
appears. The same occurs for higher T. This qualitative 
agreement with the experiment is corroborated by the evo¬ 
lution of ((V/i) 2 }, displayed in Fig. [3b- Comparing these 
results with the experimental ones (inset of Fig. [5b), one 
can confirm that the interplay of GBs barrier relaxation 
and initial conditions (initial r increasing with T), in fact, 
explains the CdTe/Si(001) local evolution. Despite this 
agreement, we remark that this model does not captures 
all aspects of the microscopic dynamics of the system as 
well as of the complex packing of polycrystalline grains. 

Coarse-grained dynamic. — Figure Eb presents the 
global roughness \W{t) = wi oc (l = L)\ versus time, which 
is expected to scale as W ~ t 3 0E]. From the linear fits 
in Fig. EJq one obtains ft = 0.51(4), 0.41(5) and 0.21(5) 
for T = 150, 200 and 300 °C, respectively. This last value 
is consistent with the KPZ one ((3k pz ^ 0.24), as was also 
found for T = 250°C pjQ. In turn, for T = 150°C, the 
value is consistent with an uncorrelated growth (where 
P = 1/2 00 ), whereas for T = 200 °C the exponent 
does not correspond to any known UC, possibly due to 
crossover effects. Indeed, for short times, one finds an 
initial slope ft e // ~ 1/2 and a tendency of /3 e // to decrease 
in time. 

From the first zero of the slope-slope correlation func¬ 
tion (SSCF) T(l,t) = (V/i(x, f)V/i(x + 1 , t)) one may esti¬ 
mate the average grain width (M) [31] . This quantity is 
expected to be of the same order of the correlation length 
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Fig. 4: (a) Global roughness W and (b) first zero of SSCF ( M ) 
versus time. The lines indicate linear fits used to extract the 
exponents {3 [in (a)] and 1/z [in (b)]. 


£(i), which scales as £ ~ t 1 / 2 , where z is the dynamic ex¬ 
ponent [B]. For all T analyzed here, T(l,t) presented an 
oscillatory behavior similar to that found in Ref. m for 
T = 250 °C. The values of ( M(t)) extracted from those 
plots are depicted in Fig. [4jo. For T = 150 °C, one finds 
£ ~ (M) const , confirming that the growth is uncorre¬ 
lated. A similar behavior is found for T = 200 °C at short 
times. In contrast, for large t, a finite z arises, namely 
z = 1-8(5), which agree with the KPZ value ( zkpz ~ 1.62) 
within the error bar. Finally, for T = 300 °C, linear fits 
at short and long times yields z « 3.0 and z ss 1.3, re¬ 
spectively. However, in this case one can not ensure that 
£ ~ (M), because small grains at the top of the multi- 
peaked mounds (see Fig. hr can make (M) smaller than 
£. This deposes against the reliability of a large z. Any¬ 
way, one notice that /3 = 0.21(5) and z ~ 3.0 are consistent 
with the Villain-Lai-Das Sarma (VLDS) class [32j, while a 
small z suggests KPZ growth. Thus, at this point, based 
only on the DSA, one can not decide what is the UC of 
the films grown at T = 200 and 300 °C. 

Beyond the scaling exponents, height distributions 
(HDs) pigi lBBHSBj . squared local roughness distributions 
(SLRDs) j36]|37]j maximal relative height distributions 
(MRHDs) [38] are also expected to be universal. In¬ 
deed, the universality of these distributions, in (2 + 1) 
KPZ class, has been experimentally demonstrated, by 
us, in the growth of CdTe/Si(001) at T = 250 °C [TT] , 
More recently, Halpin-Healy and Palasantzas have ap¬ 
plied the same method to confirm KPZ growth in oligomer 
films [T2] . We remark that SLRDs and MRHDs are built 
by measuring the squared roughness ( w 2 ) and the relative 
maximal height m = h ma x ~ (h) into boxes of lateral size 
l spanning the whole surface. The size l must be larger 
than the pixel size and l <C £ m- 

Figures [SJa-c show the HDs for all T studied. When 
T = 150 °C, experimental HDs are well-described by a 
Gaussian, as expected for an uncorrelated growth. For 
T = 200 °C, at short times, the HDs are also close to a 
Gaussian but, for large t, a nice collapse with the KPZ 
distribution is found. This agreement is confirmed by the 
skewness S = 0.43(5) and kurtosis K = 0.5(2) of the HDs, 
very close to the KPZ values S = 0.42(1) and K = 0.34(2) 



Fig. 5: Rescaled HDs for films grown at (a) T = 150, (b) 200 
and (c) 300 °C. Rescaled SLRDs for films grown at (d) T = 200 
and (e) 300 °C. (f) Rescaled MRHDs for large deposition times. 

Here, ax = \J (A 2 ) — ( X ) 2 . Insets display the same data of 
main plots in linear scale. 


SUES]. These results give further evidence of a random- 
to-KPZ crossover, a subject of wide theoretical interest 
(see [52] and references therein). Finally, for T = 300°C 
one still finds a reasonable agreement between the exper¬ 
imental HDs and the KPZ one, but now with a heavier 
left tail than the right one, yielding a negative skewness 
S = —0.2(2), with kurtosis K = 0.3(2). These values are 
consistent with KPZ class (with A < 0 in Eq. [T]) within the 
error bars. Moreover, the HDs in Fig. [5t discard VLDS 
as the possible asymptotic UC. 

The experimental SLRDs for T = 200 °C (Fig. [5J1) also 
deviate from the KPZ distribution at short times, but have 
a nice agreement for long t, giving a final confirmation of a 
crossover towards KPZ. For higher T, the SLRDs exhibit 
a nice collapse with the KPZ one (see Fig. [5^). We must 
remark that the stretched exponential decay in SLRDs 
right tail is a hallmark of the KPZ class and contrasts with 
the Gaussian decay of the VLDS distribution. Finally, the 
MRHDs for T = 200 and 300 °C are presented in Fig. [5f 
providing additional proof that CdTe grows according to 
KPZ equation. 


P-4 






















































Temperature effect on (2 + 1) experimental KPZ growth 


The relation between local and coarse-grained dynam¬ 
ics can be understood as follows. For T = 150 °C, the 
low diffusion and the energy barrier at GBs prevents co¬ 
alescence/packing of grains and, thus, the propagation of 
correlations at interface, so that inter-grain fluctuations 
evolve uncorrelated. This also happens at short times, for 
T = 200 °C, but the relaxation process at GBs gives rise 
to an asymptotic correlated growth. For higher T, where 
diffusion is more operative, these processes start early as 
well as the KPZ scaling. In terms of KPZ equation (Eq. 
[H), the random growth at low T implies v ss 0 and A ss 0. 
For T = 200 °C, one expects A > 0, but small, so that 
growth is dominated by noise initially and by non-linear 
effects asymptotically. The absence of a crossover when 
T = 250 °C [H] indicates a larger A > 0. Thus, A (T) 
seems to be a positive increasing function in this range 
of T. As discussed in Ref. nu, the possible origin of this 
KPZ growth is the complex coalescence/packing dynamics 
of the polycrystalline grains, where some grains cover their 
neighbors. Due to shape constraints, they do not necessar¬ 
ily fills all available space in its neighborhood, producing 
a positive velocity excess (A > 0) in the growth, similar 
to the lateral aggregation in ballistic deposition [5]. This 
process is more operative as higher is T, due to the larger 
inter-grain diffusion, so larger should be A. Otherwise, for 
T = 300 °C, the negative skewed HDs reveals A < 0, which 
is typical of KPZ systems where there exists deposition 
refuse as, for example, in the RSOS model [40] . Therefore, 
a possible explanation for A < 0 is that the sticking coeffi¬ 
cient is smaller in regions with very large slopes at surface. 
Indeed, one sees in Fig. [2].r that ((V/i) 2 ) for T = 300 °C is 
larger than for lower T. This can explain why this effect 
appears only at high T USE], Anyway, it is astonish¬ 
ing that so contrasting KPZ mechanisms can emerge in 
CdTe growth and indicate the possibility of control and 
even turn off the non-linearity (i.e., to make A = 0) by 
only adjusting T. 

Final remarks. — We finish stressing that the de¬ 
tailed morphological analysis performed here is impera¬ 
tive to determine the UC of the system. Since the com¬ 
plex mound evolution gives rise to finite-size corrections, 
crossover effects/transient anomalous scaling, it is not 
possible to drawn any conclusion about the asymptotic 
growth dynamic based only on the traditional study of 
the scaling exponents. This should explain why reliable 
experimental evidences of KPZ and other classes are so 
rare. Notwithstanding, we show that HDs, SLRDs and 
MRHDs are less susceptible to the above effects and have 
allowed us to determine, conclusively, that CdTe surface 
fluctuations for films grown at T £ [200, 300 °C] evolve 
according to the KPZ equation. Therefore, rather than 
a complementary analysis, the study of distributions is a 
crucial tool to unveil the growth dynamics. We believe 
that this findings will motivate future works in the same 
vein, as well as the application of these methods in previ¬ 
ously studied systems. 
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